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ABSTRACT 

The large number of exoplanets found to orbit their host stars in very close 
orbits have significantly advanced our understanding of the planetary formation 
process. It is now widely accepted that such short-period planets cannot have 
formed in situ, but rather must have migrated to their current orbits from a for- 
mation location much farther from their host star. In the late stages of planetary 
formation, once the gas in the proto-planetary disk has dissipated and migration 
has halted, gas-giants orbiting in the inner disk regions will excite planetesimals 
and planetary embryos, resulting in an increased rate of orbital crossings and 
large impacts. We present the results of dynamical simulations for planetesimal 
evolution in this later stage of planet formation. We find that a mechanism is 
revealed by which the collision-merger of planetary embryos can kick terrestrial 
planets directly into orbits extremely close to their parent stars. 

Subject headings: celestial mechanics - methods:numerical - planets and satel- 
lites:formation - stars:individual (OGLE-06-109L, 47 Ursae Majoris) 



1. INTRODUCTION 



To date over 490 extrasolar planets have been discovered, revealing a wide diversity of 



planetary systems (http://exoplanet.eu). One of more unusual phenomena so revealed has 
been the population of "Hot Jupiters" - gas-giants found in very small orbits (periods < 8d) 
about their p arent stars - of which th e prototype was the very first gas-giant exoplanet discov- 
ered, 51 Peg (IMayor fc Quelozlll995l ). It is beheved that such short-period gas-giants cannot 
have formed this close to their parent stars , and so rnust have migrated in, or been scattered 



in, from a more distant formation region (ILin et al.lll996l : IWeidenschilling &: Marzarilll996 
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Ida & Lin 


2004; 


Chambers 


2009) 



The measurement precisions that make the detection of 



such short-period exopl anets possib 



Dopp ler (e.g. Gl 876 d (IRivera et al. 



e hav e over rece nt years continua lly improv ed for both 



2005h. Gl 581 c (lUdrv et al.ll2007f) . 61 Vir b f IVogt et al. 



20101)) and transit (e.g. Kepler-4b (IBorucki et al.ll2010l )) detection. What then are the pos- 
sible formation mechanisms that can produce such close-in terrestrial and super-terrestrial 
planets? 

Several models h ave been proposed for the formation of close-in terrestrial planets. 
Raymond et al.l (120061 ) have shown that super-Earths could form interior to a migrating 
Jovian planet. As they migrate inward, such gas-giants can shepherd planetary embryos 
interior to their orbits, w hich can then further collide and merge to generate Earth-like 
planets (IZhou et al.l l2005l ). It has also been suggested that orbital mig ration and planet- 
plane t scattering could potentially produce short-period super-Earths (iBrunini fc Cionco 
20051 : iTerquem fc Papaloizoul 120071 : [Raymond et al.l |2008| ) . Whatever the mechanism for 
their formation, it is likely that such planets are common around at least low-mass stars 
jKennedv k KenvonI [2008h . 



In all these scenarios, the formation of short-period Earth-like planets is associated with 
the migration of gas-giant planets. According to the core accretion paradigm for planetary 
formation, the isolation cores in the terrestrial planet formation region, and the s olid cores 
of gas-giants, are bot h formed within ~ 1 Myr from kilometer-sized planetesimals (ISafronov 



1969 



Wetherill 



iKokubo fc Ida 



1980 



2002 



Subsequently m assive solid cores accrete disk gas to form giant planets 



Ida fc Linll2004j ) at ~ 3 — 6 Myr, before the disk disperses ( iHaisch et al. 
200 ll ). In the late stage of planet formation, when giant planets have ceased migration 
after the gas disk clears, the disk of countless planetesimals and planetary embryos will 
become turbulent due to stirring by gas-giants over hundreds of million years (or potentially 
even longer). In the meantime, it is expected that orbital crossings and g iant impacts will 



frequently occur, which could lead to t he formation of terrestrial planets (jChambersI 12001 



Raymond et al.ll2004j : IZhang fc Jill2009[ ) and short-period Earth-like planets 



In this Letter, we present a potential new formation mechanism for short-period Earth- 
like planets in the late stage of planet formation through a collision-merger scenario. In 
this mechanism, a planetary embryo is directly kicked to a close-in orbit after a collision 
with another embryo, and then the larger merged body is seized by the central star as a hot 
Earth-like planet. 
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SIMULATION SETUP 



Extrasolar planetary systems that harbor pairs of Jupiter-to-Saturn-mass companion s 



are of particular interest to researchers (jGozdziewski 112002 



e.g.. OGLE-06-1 09Lbc flGaudi et al.ll2008r) . 47Umabc flButler fc Marcvlll996 : Fischer et al. 



Ji et al 



2005 



Zhang et al. 2010) 



20021 ). G1876 be (IMarcy et al.ll200l[ ). It is interesting to consider whether it is likely that such 
systems might host addi tional hot terrestrial planets (as, for example, the G1876 system does 



in the form of G1876 d - iRivera et al.l ( l2005l )). and further how such planets might form and 



evolve. We have therefore carried out simulations that explore such a system architecture. 

In t otal, 30 runs were performed using a hybrid symplectic algorithm in the MERCURY 
package (jChambersI Il999l ) for following two such systems. The initial conditions of the two 
systems simulated were: 



Simulation 1 - two giant planets are simulated with initial orbital parameters (Mp, a, e^) 
= (0.71 Mjup, 2.3 A U, 0.001) and (0.2 7M.i„p, 4.6 AU, 0.11), to emulate the OGLE-2006- 
BLG-109L system ( Gaudi et al. 2008). 500 planetary embryos and planetesimals0 with 
total mass 10 were distributed between 0.3 AU < a < 5.2 AU and with e < 0.02. 
Each of the 26 runs carried out over 400 Myr. 

Simulation 2 - two giant planets are simulated with initial orbital parameters (Mp, a, Cp) 
= (2.9 M.T„p, 2.08 AU 0.05) and (l.lMj^p, 3.97 AU, 0.001), to emulate the 47Uma 
system (IFischer et al.ll2002l ). 648 planetary objects with total mass of 5.14M®were 
distributed in the region 0.3 AU < a < 1.6 AU with e < 0.02. Each of the four runs 
evolved over 100 Myr. 



The other initial orbital elements of each planetary embryo (or planetesimal) are ran- 
domly generated - the arguments of periastron, longitudes of the ascending node, and mean 
anomalies range from 0° to 360°, and inclinations are from 0° to 1°. In addition, the hybrid 
integrator parameters are adopted as a stepsize of 3 days (~ a twentieth of a period for the 
innermost possible body at 0.3 AU), and a Bulirsch-Stoer tolerance of 10~^^. At the end of 
integration, the changes of energy and angular momenta are 10~^ and 10"^^, respectively. 
In these runs, the gravitational interactions of all bodies are taken into account. Two bodies 
are assumed to collide whenever the distance between them is less than the sum of their 



physical radii (jChambersI Il999[ ) . If two objects collide, they are merged into a single body, 
without fragmentation, after the collision. 



^Herein the masses of embryos range from several lunar- mass to Mar-mass, and those of smaller "plan- 
etesimals" have approximately a lunar mass, rather than a planetesimal mass. 
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3. RESULTS 

3.1. Simulation results 

In our simulations, we find that the collision-merger mechanism produces close-in ter- 
restrial planets in 20% of the runs carried out (5 of 26 Simulation 1 runs, and 1 of 4 
Simulation 2 runs). Th e simulations exhibit a classical plan etary accretion scenario in 



their late stage formation (jChambersll200ll : [Raymond et al.ll2004j ). Figure [H shows snapshots 



at various evolution times for a representative run of Simulation 1. Initially, the embryos 
and planetesimals reside in a cold disk, which is quickly stirred by the two gas-giants and 
excited to highly eccentric orbits within 0.1 Myr. We also see that three small bodies are 
involved in a 1:1 resonance with the inner giant by that time. By the end of IMyr, most of 
the initial objects have been removed by ejection or collision due to frequent orbital crossings 
in this chaotic stage. In addition, we see that a close-in planet has formed at ~ 1 Myr which 
subsequently remains very stable. At the conclusion of the run (400 Myr), three planetes- 
imals survive, of which one has been seized as a Trojan body by the inner giant, and the 
other two move at ~ 1 AU in eccentric orbits. 

Figure [2] shows the time evolution of the mass, semi-major axis, and eccentricity of 
the short-period terrestrial planet formed in the Simulation 1 run shown in Fig. [1] At 
0.0356 Myr, two bodies that may be excited by secular resonance of gas-giants, collide at 
very high eccentricities (e = 0.91 and 0.80, shown by the red and black lines in Fig. |2l 
respectively) and are then assumed to merge into a single planetary embryo. That merged 
body (the remaining black line in Fig. is captured by the parent star as a short-period 
planet, and its orbit dramatically shrinks from ~ 0.4 AU at the time of the collision, down 
to 0.077 AU. Subsequently, three additional collisions take place over the further late-stage 
evolution of that merged object. Fig. [2] shows that the embryo moves slightly inward at 
each collision, and that its mass also increases. Moreover, we note that it finally becomes 
a 3.3 Mercury-mass planet with a close-in orbit about 0.056 AU, and its eccentricity drops 
down to e=0.13 after the last collision. The orbit may then, of course, be further circularized 
by tidal interaction with the star over even longer timescales. 

Figure [3] shows the formation and evolution of a similar terrestrial planet that emerges in 
one of the runs for Simulation 2. At 2.2 Myr, the semi-major axis of a ~ 0.9 Mercury- mass 
embryo drops down from ~ 0.8 AU to 0.06 AU as a result of a collision with a highly-eccentric 
planetesimal excited more than a million years earlier. The merged body has an eccentricity 
that drops from 0.90 to 0.50 immediately after the impact. The enlarged, merged body 
subsequently undergoes additional collisions, and its eccentricity further evolves to e=0.33 
(after its last collision) with a final mass of 1.3 Mercury-mass. Here, the collision- merger 
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Fig. 1. — A snapshot of planet formation in the late stage for Simulation 1. The panels show 
the orbital eccentricity versus semi-major axis for each surviving body at simulation times 
of 0, 0.1, 1.0, 10, 50 & 400 Myr. The radii and the color of the embryos and planetesimals 
are related to their mass, with radius proportional to m^^^. The two giants are, respectively, 
at 2.3 and 4.6 AU. A close-in terrestrial planet forms at ~ IMyr and it remains stable over 
secular evolution. 
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Fig. 2. — Mass, semi-major axis and eccentricity evolution of the short-period terrestrial 
planet the emerges from the Simulation 1 run shown in Fig. [H The black and red lines in 
the lower two panels show the semi-major axis and eccentricity evolution of the two bodies 
that collide to form a merged planetary embryo, which is kicked from 0.4 AU to 0.077 AU at 
0.0356 Myr. Subsequently that merged embryo (shown as a single black line after 0.0356 Myr) 
is subject to further collision-mergers, with the epoch of each collision shown by the solid 
triangles. The resultant mass evolution of this body is shown in the upper panel. 
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Fig. 3. — Mass, semi-major axis and eccentricity evolution of a short-period terrestrial planet 
that emerges from a run of Simulation 2 - layout is the same as for Fig. O In this case the 
first collision-merger event occurs at 2.2 Myr, and the embryo is thrown from its location of 
~ 0.8 AU at the time of the collision to 0.06 AU. 
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scenario may provide some clues of the origi ns of the moderate eccentricities seen in super- 
Earths detected to date (e.g., HD 181433 b feouchv et aPboOoh l 



The major difference in the evolution of these two examples is that the short-period 
planet that evolved in Simulation 1 was moved into an inner orbit at a very early stage, 
and subsequently accreted a majority of the mass available in nearby orbits; while the 
Simulation 2 planet had almost completed accretion into a Mercury-mass embryo before 
it moved closer to the star. In all simulations, we notice that terrestrial planets and bodies 
formed at short periods via collision-merger events come into bei ng within 10-30 Myr, which 
agrees with the estimated timescale of terrestrial core formation (lYin et al.ll2002l ). as derived 
from the chronornetry of meteorites and numerical sirnulatio ns of terrestrial planet formation 
dChambersI boOll : iRavmond etaP bood : Izhang fc Jil boooh . In addition, we find that all 
survivors remain stable in their final configurations. 

These results indicate that a collision-merger mechanism could indeed produce short- 
period, terrestrial planets in two systems that host two gas-giants. Similarly, we also find 
the above outcomes in other 4 runs. However, a natural question then arises - do the bodies 
that take part in these collisions really merge? Or will they become fragmented? 



3.2. Merger versus fragmentation 



In the accretion model of MERCURY, a collision-merger scenario oc curs whenever th e 
distance between two bodies is less than the sum of their physical radii ( IChambersI Il999l ) . 
and the package models the two bodies merging inelastically to form a single new body 
that conserves mass and total momentum. The collisions in the runs, therefore, are con- 
sidered to be perfect gravitational aggregations, which assumes that enough energy is dis- 
sipated in the collision for the two bodies to remain gravitationally bound. However, ac- 
tual collisions could have a result that ranges anywhere from this result (complete merger), 
through partial fragmentation, to t he complete shattering and disintegration of both im- 
pactors. ( IWetherill fc Steward Il993[ ). 



Whether these bodies either fragment or cohere in a collision will obviously depend on 
the - currently poorly understood - physical properties of the colliding bodies ( IWetherill 



19801 ). What can be said is that the outcome will be extremely complicated. To assess the 
likely state of the merger vs fragmentation for two bodies in a collision, we can, though, 
make order-of-magnitude estimates. 



Consider two bodies of the same mass m, with relative velocity at infinity cr, and the 
sum of the physical radii Rg. The collision velocity for a head-on collision between them is 
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flSafronovlll969l : IWetheril]lll980l : lArmitagel 120071 ). 



where Vesc= \/ 4:Gm / Rs is the escape velocity at the point of colhsion, a parameter used to 
evaluate whether they will physically collide. Take the coefficient of restitution as e, then 
accretion will result if eVc < Vgsc, even if the initial impact results in fragmentation into two 
bodies. Conversely, the bodies will be unbound if 6f ^ > Vesr . - Thu s, the threshold value of 



the coefficient of restitution for these outcomes is (lArmitagd 120071 ) 

-1/2 



1 + 



a 



(2) 



This shows that if a <^ Vesc, merger and growth is likely unless collision is totally elastic; 
whereas a 3> Vesc leads to fragmentation. 

For the Simulation 1 run shown in Fig. |2] , we have used the above equations to 
assess the outcome of the ffist collision as it happens between two identical Mercury-like 
embryos, which allows us to make a rough evaluation of the likely outcome by calculating 
the instantaneous velocities of the impactors at the epoch just before the collision. Now we 
notice that at the ffist collision the body was impacted onto a close-in orbit. The Vesc of the 
two impactors are nearly the same - 3.12 kms~^ (assuming equal bulk density); the velocities 
of the impactors at the collision epoch near the pericenter are estimated to be 43.21 km s~^ 
and 35.50 km s~^, respectively, thus we have an approximate relative velocity projected to 
the relative position of two colliding bodies of 13.49 km s~^. In this case, a merger requires 
e < 0.23, which is close to the accretion co ndition (e < 0.34) in realistic accretion model for 
head-on collision ( iKokubo &: Gendal 120101 ). On the basis of above analysis, a merger seem 
to be possible for two eccentric objects when the collision occurs in the nearby region of 
central star, subsequently the merged body is seized by the star at close-in orbit. In the 
collision-merger process, moderate energy should be released, and they could be converted 
into the internal heat of the merger in the collision between embryos, e.g., simulations of a 
supposed Moon-forming impact show that the collision can deliver pr odigious ener gy to the 
Earth, which could lead the proto-Earth to a mixed solid-melt state ( ICanupl 120081 ) Ji 

We also obtain similar resu lts for the Simulation 2 run shown in Fig. [31 In addition, 
Leinhardt fc Richardson! ( 120021 ) showed that a large mass ratio between two impactors will 
tend to lead to merger and aggregation - the accretion probability is ~ 60% (averaged over 



^At the very time before/after the coUision, the fractional energy change due to integrator was about 9 
part in 10^. Additional en ergy loss may arise from the ejection of other embryos or transfer to the envelope 



and core of giant planets ( Li et al. 20ld} 



- 10 - 



all impact parameters) for average mass ratio of 1:5. This suggests that the first collision 
seen in this run, where the mass ratio of 1:3.43, is likely to result in a merger. 



4. DISCUSSION and CONCLUSION 

We have uncovered a new mechanism for producing short-period terrestrial planets via 
collisions-mergers in the late stages of planetary formation. In this mechanism, two highly- 
eccentric bodies first undergo a severe orbital crossing and then form a short-period planet 
via collision-merger. In the set of simulations performed to date, this mechanism produces 
a short-period, terrestrial planet in 20% of runs. 

As mentioned previously, the formation rate for short-period terrestrial planets via a 
collision-merger process is only a moderate 20%. However, this low rate may be a result of the 
limits imposed on our simulations by current computational capabilities, which restrict our 
adopted population of embryos and planetesimals to a few hundred objects with a total mass 
of only several times that of the Earth. The resultant planetesimal disk in our simulations 
is much smaller than that of the Minimum Mass Solar Nebula (~ 0.01 Mq within 30 AU 



( iWeidenschillingl 119771 : lHayashilll98ll )) - which would also contain billions of small bodies. 
Increasing the number of bodies and the total mass of the proto-planetary disk would likely 
increase the efficiency with which this mechanism produces short-period terrestrial planets. 

In addition, it is worth noting that close-in planets emerge from our simulations within 
a few million years. This is a significantly shorter timescale than the billion years over 
which the Solar System is thought to have undergone significant evolution. So, while near- 
infrare d observations of y oung cluster samples, indicate an overall dust disk lifetime of ~ 



6Myr (iHaisch et al.ll200ll ). the planetary system will actually continue to evolve over much 
longer timescales following the clearing of the dust and gas disk. During this late stage 
of planetary formation, frequent orbital crossings and huge impacts will occur, which are 
likely to significantly boost the feasibility of collision-merger events producing short-period 
terrestrial bodies. 

The collision-merger scenario for the formation of short-period planets does not require 
perfect accretion. Rather it relies on the collisions pushing the resultant body inward, so 
that the central star can capture it as a short-period planet. In this sense, such a mechanism 
could play a key role in throwing the largest fragments resulting from severe impacts into 
short-period orbits. On the other hand, given the diversity in the architectures of currently 
known systems, exoplanets are likely to form through a variety of mechanisms rather than 
through a uniform dominant process (D. Lin 2009, private communication). Our simulation 
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results show one potential mechanism for the origin of short-period terrestrial planets in a 
compact disk with two gas-giants, and may predict an abundance of close-in bodies for this 
family. 
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